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ABSTRACT 

We investigate the scaling relations between the supermassive black holes (SMBHs) mass 
(M bb ) and the host bulge mass in elliptical galaxies, classical bulges, and pseudo-bulges. We 
use two-dimensional image analysis software BUDDA to obtain the structural parameters of 
57 galaxies with dynamical M bb measurement, and determine the bulge X-band luminosi- 
ties (Xm.kX stellar masses (M s ), and dynamical masses (M d ). The updated M bh -L bulK , M bh - 
Af s , and M bh -M d correlations for elliptical galaxies and classical bulges give A/ bh ~0.006A/ s 
or 0.003 M t . The most tight relationship is log(M bh /M ) = a + /31og(M d /lO n M ), with 
a = 8.46 ± 0.05, (3 — 0.90 ± 0.06, and intrinsic scatter eo = 0.27 dex. The pseudo-bulges 
follow their own relations, they harbor an order of magnitude smaller black holes than those 
in the same massive classical bulges, i.e. Af bb ~0.0003A/ s or 0.0002Af d . Besides the M bh -cr* 
(bulge stellar velocity dispersion) relation, these bulge type dependent Af bh -M bul scaling rela- 
tions provide information for the growth and coevolution histories of SMBHs and their host 
bulges. We also find the core elliptical galaxies obey the same Af bh -M d relation with other 
normal elliptical galaxies, that is expected in the dissipationless merger scenario. 

Key words: black hole physics - galaxies: bulges - galaxies: formation - galaxies: funda- 
mental parameters - galaxies: nuclei. 



1 INTRODUCTION 

Supermassive black holes (SMBHs) are now believed to be a key 
element in galaxy formation. They grow and coevolve with their 
host galaxies, regulate star formation and heat intracluster medium. 
The most indicative evidences of this symbiosis are the tight corre- 
lations between the SMBH masses (M bh ) and bulge properties. In 
their early review, Kormendy and Richstone (1995) found the Af bb - 
L bld (bulge luminosities) and the equivalent M bh -M s (bulge stellar 
masses) correlation for a sample of eight SMBHs in the local qui- 
escent galaxies. For an enlarged sample of 32 galaxies, Magorrian 
et al. (1998) confirmed the M bb -L bu i and found M bh -AI d (bulge dy- 
namical masses) correlation with M bb ~0.006Af d and an intrinsic 
rms scatter eo < 0.5 dex. The M bb -M bu i relations are often called 
"Magorrian relations" in the subsequent literature. Similar correla- 
tions are found in active galaxies (e.g., Wandel 1999, 2002; McLure 
& Dunlop 2001; Bentz et al. 2009; Gaskell & Kormendy 2009). 

The two most widely used versions of the Magorrian rela- 
tions in local quiescent galaxies are determined by Marconi & Hunt 
(2003, hereafter MH03), and by Haring & Rix (2004, hereafter 
HR04). To avoid the effect of large variations of mass-to-light ra- 
tio and intrinsic dust extinction of bulges in optical (i.e. B or K) 
bands in the previous studies, MH03 used near-infrared image to 



E-mail: jhu@mpa-garching.mpg.de 



measure the bulge properties. They obtain the structural parameters 
of 37 galaxies by a two-dimensional bulge/disk decomposition pro- 
gram GALFIT (Peng et al. 2002), determined L bu] and estimated Af d . 
They found tight M bh -L bu i and M hh -M L] relations (M bh ~0.002A/ d , 
eo — 0.3 dex). In order to measure the bulge masses more accu- 
rately, HR04 derived M d through solving Jeans equation of dynam- 
ical model for a sample of 30 galaxies (including 12 from the litera- 
ture), and determined a similar tight M bb -M d relation (Af bh -~0.0014 
M it e ~ 0.3 dex). 

Although the M bb -Af bl ,i correlations have been well defined as 
ones of the most tight SMBH-bulge scaling relations (Novak et al. 
2006), several important problems are still unclear. 

In a previous study, we have demonstrated the A/ bh -<7» (central 
stellar velocity dispersion) relations are different for classical 
bulges and pseudo-bulges (Hu 2008, hereafter H08). The pseudo- 
bulges are a kind of disk-like bulges in the center of disk galaxies, 
they have distinct properties and origins from the classical bulges 
(e.g., Kormendy & Kennicutt 2004). On average, the SMBHs in 
pseudo-bulges are ~6 times smaller than in their classical counter- 
parts. Obviously, it is worth exploring whether the similar differ- 
ences exist in the Af bb -M blI i relations. 

Kormendy & Gebhardt (2001) once checked the B-band M bh - 
L bu] relation for five pseudo-bulges, and found on significant in- 
consistency with that of the classical bulges. However, their re- 
sult was not conclusive due to the small sample, the wrong bulge 
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type identification of a galaxy (NGC 4258), and the large intrin- 
sic scatter of the B-band M bh -L bu] relation. Recently, Greene et al. 
(2008) found the M bb /M bu , ratio of a sample of active black holes 
in pseudo-bulges and spheroidal galaxies is about an order of mag- 
nitude lower than that in local quiescent classical bulges. However, 
it is not clear whether the difference come from the activities of 
SMBHs or redshift evolution (most of their sample galaxies are of 
redshift z = 0.1 ~ 0.2). Assuming two kinds of bulges follow the 
same M bh -Af s relation, Gadotti & Kauffmann (2009) found pseudo- 
bulges have different Af bb -<7* relation, confirming the result of H08, 
but their assumption should be checked. 

Some elliptical galaxies have central core shape surface 
brightness profiles (e.g., Lauer et al. 1995). The cores are believed 
to be formed due to the interaction of binary SMBHs with sur- 
rounding stars in the dissipationless (dry stellar) mergers events. 
The brightest cluster galaxies (BCGs) and brightest group galaxies 
(BGGs) are the largest ones of core ellipticals. In H08, we have 
found the Mw(f* relation of the core elliptical galaxies are slightly 
steeper than that of the normal elliptical galaxies. On the other side, 
Lauer et al. (2007a, hereafter L07a) suggest the Ad bh -L bal relations 
are the same for all ellipticals galaxies, which deserve further ex- 
amination. 

The M bh -L bul (M d ) relations have been updated since the work 
of MH03 and HR04. Graham (2007, hereafter G07) corrected sev- 
eral issues in the previous work, adjusted the B, R and A-band 
M bb -L M relations. Gultekin et al. (2009b, hereafter G09b) im- 
proved the l/-band M bh -Z/ blll relation based on a sample of 49 galax- 
ies. Although the results of different authors derived by different 
samples are roughly consistent with each other, pseudo-bulges and 
classical bulges are mixed to fit the relations in all the previous 
work. As we have shown in H08, the slopes (/3) and intercepts (a) 
of black hole-bulge correlations are very sensitive to the low mass 
(e.g., pseudo-bulges) and high mass (e.g., core elliptical galaxies) 
objects in the sample. If the M bb -M blll relations for pseudo-bulges 
are below that for classical bulges (just like the Af bb -cr* relation), 
then a will be underestimated and /3 be overestimated. 

The local SMBH mass function (BHMF) can be derived by 
combining the black hole-bulge scaling relations (i.e. M bb -a, , M bb - 
Af bll] ) and galaxy <7„, M M distribution functions (e.g., Shankaret al. 
2004, 2009; Marconi et al. 2004). In these calculations, the shape 
and normalization of BHMF (and the local SMBH density) are re- 
lated to a and /3, the high and low end of BHMF are very sen- 
sitive to the intrinsic scatter of these relations. Another important 
parameter, the average radiative efficiency of black hole accretion 
in AGNs derived by the Soltan (1982) argument, is also sensitive to 
a and f3 (e.g., Yu & Tremaine 2002; Marconi et al. 2004; Shankar 
et al. 2004; Yu & Lu 2008). 

Considering the above mentioned observational and theo- 
retical importance, the previously determined A/ bb -M bu i relations 
should be re-examined. Therefore, in this paper, we investigate 
these correlations for different bulge types in an up-to-date sample. 
The paper is organized as follows. We describe the sample and data 
analysis processes in section 2. The black hole-bulge properties re- 
lations are shown in section 3. Finally we summarize and discuss 
our results. 

Throughout this paper, we use the base 10 logarithms, and the 
cosmological parameters Q M = 0.3, £Ia = 0.7, Ho = 70 km s^ 1 
Mpc- 1 . 



2 DATA 

2.1 Sample 

We select a sample of galaxies with secure A/ bh measurements (by 
stellar or gas dynamics, masers or stellar orbit motions). The sam- 
ple is based on that used in H08, with the following updates. 

Cygnus A (Tadhunter et al. 2003), Fornax A (=NGC 1316) 
(Nowak et al. 2008), Abell 1836-BCG (=PGC 49940), Abell 3565- 
BCG (=IC 4296) (Dalla Bonta et al. 2009), NGC 3585, NGC 3607, 
NGC 4026, NGC 5576 (Gultekin et al. 2009a), NGC 524, NGC 
2549 (Krajnovic et al. 2009) are added to the sample. The A/ bb and 
a, of these objects are taken from the literature. The la error of 
M bb for NGC 524 and NGC 2549 are estimated by the author from 
the Figure 8 of Krajnovic et al. (2009). 

M bh of NGC 821, NGC 3377, NGC 3379, NGC 3384, NGC 
3608, NGC 4291, NGC 4473, NGC 4564, NGC 4649, NGC 4697, 
NGC 5845, NGC 7457 taken from Gebhardt et al. (2003) are in- 
creased by 9% due to an numerical error in the original published 
version (G09b). 

As in H08, the distance to galaxies are taken from the mea- 
surement of surface brightness fluctuation (SBF) method, e.g., 
NGC 1316, IC 4296 (Jensen et al. 2003), NGC 3585, NGC 3607, 
NGC 4026, NGC 5576, NGC 524, NGC 2549 (Tonry et al. 2001); 
or Hubble recession velocities corrected for Virgo centric infall, 
e.g., Cygnus A, and PGC 49940. Their M bh given by the literature 
are modified accordingly. 

According to the SINFONI observation on the central stellar 
kinematics of NGC 3227 (Davies et al. 2006), its bulge locates be- 
low the oblate rotator line in the (V m /cr*)-e diagram (e.g., Binney 
1978), where V m is the maximum line-of-sight rotational velocity 
of the bulge, e is the ellipticity of the bulge. Our fit of the Sersic 
index n — 2.0 (see below) is also much larger than n = 1.1 given 
by Gadotti (2008). Therefore, we treat it as a classical bulge in this 
paper. The ct„ of NGC 3227 is modified to 131 km s _1 (Onken et 
al. 2004). 

As mentioned in H08, NGC 2787 and NGC 3384 contain both 
pseudo-bulges and classical bulges (Erwin 2008). In the following 
analysis, we treat these two objects neither as classical bulges nor 
as pseudo-bulges, but compare them with other objects. 

We include the milky way as a pseudo-bulge in the following 
analysis. The parameters of the milky way are M bb = (4.1 ± 0.6) x 
10 6 M Q (Ghez et al. 2008), bulge Sersic index n = 1.3, effective 
radius R e = 0.7kpc,a»= 103±20kms _1 (Tremaine et al. 2002), 
A'-band luminosity of the bulge log(L bu i_ K /L0 iK )= 10.25 ± 0.30 
(Dwek et al. 1995; MH03), A'-band luminosity of the total galaxy 
is log(I/ bu i, K /L 0iK )= 10.94 (Drimmel & Spergel 2001), bulge stellar 
mass M s = (1.3 ± 0.5) x 10 10 M (Dwek et al. 1995). 

Our sample consists of 58 galaxies (including the milky way), 
28 are elliptical galaxies, 22 are disk galaxies with classical bulges, 
6 are pseudo-bulges, 2 have both classical bulges and pseudo- 
bulges. 

2.2 Image decomposition 

Like MH03, we use the A-band images of samplegalaxies from 
the Two Micron All Sky Survey (2MASS) databasajfor structural 
analysis, in order to minimize the effect of mass-to-light ratio vari- 
ation and intrinsic dust extinction, especially in the pseudo-bulges 
with young stellar population and dusty structures. The images 

1 http://www.ipac.caltech.edu/2mass/ 
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have been sky background subtracted and photometric calibrated 
by the 2MASS team. The pixel size of the 2MASS image is 1.0 
arc sec. 

We use a two-dimensional bulge/disk decomposition program 
BUDD^v2.2 to measure the structural parameters of disks, bulges, 
bars and central compact sources in the galaxies. Decompostion by 
BUDDA has been tested to be very robust for components larger 
than the image PSF (e.g., Gadotti 2008, 2009). The typical PSF 
size (3 arcsec) of 2MASS images is small enough for the nearby 
galaxies in our sample, most of which are of size of several hundred 
arc sees. 

In BUDDA, the surface brightness profile of disk component is 
described by an exponential function: 

A» d (r) = mo + 1.086r/h, (1) 

where r is the galactocentric distance, fio is the disk central surface 
brightness, h is the disk scalelength. 

The surface brightness profiles of bulge and bar components 
are described by Sersic (1968) function: 

Jib(r) =^ + c n [{r/r c ) l/n - 1], (2) 

where /i c is the effective surface brightness, r c is the effective ra- 
dius, n is the Sersic index, and c n = 2.5(0.868n — 0.142) is a 
constant. 

The shape of disks, bulges and bars are modeled by concentric 
generalized ellipses (Athanassoula et al. 1990): 

(\x\/ a y + {\ y \/by = i, o) 

where x and y are pixel coordinates, a and b are the semimajor 
and semiminor axes length, c is the shape parameter. In the case 
of ellipse, c = 2, which is taken as the fixed value for disks and 
bulges. The code can fit the image with position angles, ellipticities 
(e = 1 — b/a), and c (only for bar component). If it is necessary, 
the truncation radius of the components can be fitted by the code. 

AGNs and other central unresolved sources (e.g., nuclear star 
clusters) are modeled as a point source convolved with the PSF 
Moffat profile. The accurate PSF size for each image is determined 
by fitting the foreground stars. 

The best model is achieved by comparing the % 2 of fitting with 
several sets of initial values. Finally, the code generate the model 
image and the residual image for evaluating the quality of fitting. 

The best-fit structural parameters of the galaxies are listed in 
Table 1 . The photometric parameters in Table 1 have been corrected 
with the Galactic and internal extinction. The A'-band Galactic ex- 
tinction (Ak) is taken from the NED databas^] which is calcu- 
lated based on Schlegel et al. (1998). The internal extinction due 
to the inclination of the galaxy is taken from the HYPERLEDA 
databas^] For NGC 224 and Circinus, Ak estimated by Schlegel 
et al.'s routine is unreliable, we follow the correction used by HY- 
PERLEDA. Some objects (e.g., NGC 221) are labeled as ellipti- 
cal galaxies by NED, but they are fitted better adding disk compo- 
nents. These model dependence will induce systematic difference 
between our results and the previous work (see discussion below). 

The typical la relative errors for individual structural param- 
eters (column 2-11 in Table 1) are 10%-20%, but some parameters 

2 BUDDA is developed by R. E. de Souza, D. A. Gadotti, and S. dos An- 
jos (2004). The documents and the program of BUDDA are available at 
http://www.mpa-garching.mpg.de/~dimitri/budda.html 
° http://nedwww.ipac.caltech.edu/ 
4 http://leda.univ-lyonl.fr 
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Figure 1. Four model components of NGC 1068. 

are strongly coupled (e.g., rib and r c ). The photometric uncertain- 
ties of m ga i is about 0.05 mag, the la error of the luminosity frac- 
tion of each component (column 14-17 in Table 1) is less than 5%. 
As we are more interested in the statistical properties of the sample, 
the errors of the parameters are not listed in Table 1 . 

The BUDDA can generate the model images for each compo- 
nent of the galaxy. As an example, Figure 1 presents the four com- 
ponents of NGC 1068. 

The images of our best-fit models are shown in Figure 2. Most 
of the models are visually very good to resemble the observation. 
In the residual images, structures not included in the model are 
prominent, such as spiral arms (NGC 4258), dust lane (NGC 5 128), 
optical jet (NGC 4486), and asymmetric central structure (NGC 
4486A). In some disk galaxies with large inclination (e.g., NGC 
3079), the bar component in the model may be in fact the spiral 
structures in the disc. Three edge-on galaxies (NGC 2549, NGC 
3115, NGC 4026) are not well fitted with BUDDA, we remove them 
from the sample in the following analysis. 

2.3 Bulge mass 

The bulge properties of sample galaxies are calculated based on the 
two-dimensional decomposition, the results are listed in Table 2. 

The bulge luminosity L bll]K is the the product of the total 
galaxy luminosity L M , K and the bulge luminosity fraction (B/T). 
The la error of L bu i, K is adopted as 10%. The bulge stellar mass M s 
is the product of L bulK and A"-band mass-to-light ratio M/L. We 
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Table 1. Galaxy structural parameters of K band image 



Galaxy 


/'() 


h 






r e 


n b 




Mbar 


'^bar 


^bar 


m c 




B/T 


D/T 


Bar/T 


C/T 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


N221 


14.5 


24 


0.20 


14.1 


10.5 


4.0 


0.23 


1 


/ 


1 


/ 


5.07 


0.453 


0.547 








N224 


15.2 


345 


0.55 


15.5 


118.4 


1.9 


0.14 


1 


/ 


1 


11.69 


0.95 


0.323 


0.676 





0.001 


N524 


/ 


/ 


/ 


17.8 


35.0 


4.0 


0.04 


1 


/ 




' 


7.09 


1 











N821 


17.3 


26 


0.19 


15.7 


7.9 


3.3 


0.37 


1 


/ 


1 


1 


7.87 


0.631 


0.369 








N1023 


15.8 


45 


0.63 


14.7 


7.4 


2.5 


0.17 


1 


/ 


i 


1 


6.23 


0.347 


0.653 








N1316 


15.8 


42 


0.26 


15.1 


13.3 


3.1 


0.33 


1 


/ 


1 


1 


5.57 


0.356 


0.644 








N1399 


/ 


/ 


/ 


16.7 


27.3 


3.4 


0.09 


1 


/ 


1 


1 


6.58 


1 











N2549 


17.2 


36 


0.76 


15.8 


9.7 


2.9 


0.45 


1 


/ 


1 


1 


8.14 


0.830 


0.170 








N2974 


/ 


/ 


/ 


16.9 


19.8 


3.9 


0.36 


1 


/ 


I 


1 


7.63 


1 











N3031 


15.6 


102 


0.43 


15.3 


26.0 


3.4 


0.29 


1 


/ 


i 


1 


3.78 


0.248 


0.752 








N3115 


14.7 


35 


0.70 


13.9 


6.7 


2.5 


0.43 


1 


/ 


i 


1 


5.88 


0.294 


0.644 


0.062 





N3227 


16.7 


37 


0.56 


15.7 


4.8 


2.0 


0.37 


1 


/ 


1 


13.0 


7.67 


0.138 


0.768 





0.086 


N3245 


17.3 


31 


0.55 


15.6 


8.0 


3.9 


0.36 


1 


/ 


1 


' 


7.85 


0.641 


0.359 








N3377 


16.5 


26 


0.50 


15.3 


7.7 


4.4 


0.44 


1 


/ 


1 


1 


7.43 


0.523 


0.477 








N3379 


/ 


/ 


/ 


17.2 


35.1 


4.8 


0.11 


1 


/ 


i 


1 


6.25 


1 











N3414 


17.6 


32 


0.26 


15.6 


6.8 


3.0 


0.25 


1 


/ 


1 


1 


7.94 


0.557 


0.443 








N3585 


18.2 


66 


0.45 


16.3 


20.5 


4.0 


0.44 


1 


/ 


i 


1 


6.68 


0.748 


0.252 








N3607 


/ 


/ 


/ 


17.0 


24.9 


2.9 


0.20 


1 


/ 


1 


U.2 


7.01 


0.922 








0.078 


N3608 


/ 


/ 


/ 


17.5 


18.6 


4.6 


0.18 


1 


/ 


i 


1 


8.10 


1 











N3998 


18.7 


97 


0.23 


15.6 


9.5 


4.2 


0.17 


1 


/ 


1 


1 


7.34 


0.851 


0.149 








N4026 


16.4 


42 


0.89 


16.3 


13.0 


4.7 


0.35 


1 


/ 


I 


1 


7.57 


0.765 


0.235 








N4151 


18.0 


59 


0.30 


16.5 


9.5 


4.0 


0.12 


1 


/ 


1 


11.7 


7.38 


0.403 


0.321 





0.206 


N4258 


16.2 


80 


0.50 


15.6 


10.6 


2.6 


0.37 


1 


/ 


1 


1 


5.36 


0.107 


0.893 








N4261 


/ 


/ 


/ 


17.6 


29.5 


4.1 


0.21 


1 


/ 


1 


1 


7.40 


1 











N4291 


/ 


/ 


/ 


16.4 


10.0 


4.4 


0.21 


1 


/ 


i 


1 


8.46 


1 











N4459 


16.5 


27 


0.21 


15.2 


6.8 


2.5 


0.12 


1 


/ 


1 


1 


7.10 


0.431 


0.569 








N4473 


/ 


/ 


/ 


16.4 


19.9 


4.1 


0.41 


1 


/ 


1 


1 


7.15 


1 











N4486 


/ 


/ 


/ 


17.5 


55.8 


3.0 


0.07 


1 


/ 


1 


1 


5.83 


1 











N4486A 


/ 


/ 


/ 


15.7 


6.1 


4.2 


0.32 


1 


/ 


1 


1 


8.98 


1 











N4552 


1 


/ 


/ 


16.7 


21.9 


4.6 


0.06 


1 


/ 


1 


1 


6.72 


1 











N4564 


16.2 


20 


0.65 


15.6 


7.2 


3.6 


0.38 


1 


/ 


i 


1 


7.82 


0.497 


0.503 








N4596 


17.5 


49 


0.39 


16.1 


8.0 


3.3 


0.13 


18.6 


46 


1.0 


1 


7.10 


0.277 


0.603 


0.120 





N4621 


16.2 


32 


0.45 


15.3 


8.5 


3.5 


0.32 


1 


/ 


i 


1 


6.73 


0.399 


0.601 








N4649 


/ 


/ 


/ 


17.1 


47.1 


3.4 


0.16 


1 


/ 


1 


1 


5.78 


1 











N4697 


16.2 


36 


0.36 


15.8 


13.7 


3.0 


0.43 


1 


/ 


1 


1 


6.47 


0.393 


0.607 








N4742 


18.5 


49 


0.14 


14.9 


4.6 


3.7 


0.35 


1 


/ 


1 


1 


8.33 


0.764 


0.236 








N5077 


/ 


/ 


/ 


16.5 


12.1 


3.2 


0.24 


1 


/ 


1 


1 


8.20 


1 











N5128 


15.3 


67 


0.07 


15.7 


28.9 


3.0 


0.06 


1 


/ 


1 


1 


3.87 


0.287 


0.713 








N5252 


17.1 


25 


0.30 


15.5 


6.0 


3.2 


0.21 


1 


/ 


i 


1 


7.93 


0.509 


0.491 








N5576 


/ 


/ 


/ 


17.6 


11.3 


4.0 


0.32 


1 


/ 


1 


1 


9.25 


1 











N5813 


17.4 


43 


0.15 


17.0 


12.7 


4.6 


0.12 


1 


/ 


i 


1 


7.41 


0.479 


0.521 








N5845 


/ 


/ 


/ 


15.2 


4.9 


4.6 


0.37 


1 


/ 


1 


1 


9.05 


1 











N5846 


/ 


/ 


/ 


17.8 


36.6 


3.1 


0.04 


1 


/ 


i 


1 


6.99 


1 











N6251 


/ 


/ 


/ 


17.3 


10.8 


3.4 


0.11 


1 


/ 


1 


1 


8.93 


1 











N7052 


16.8 


20 


0.59 


16.8 


9.1 


3.4 


0.29 


1 


/ 


i 


1 


8.56 
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N3079 


16.2 


47 
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0.2 


14.5 


7.01 


0.233 


0.607 


0.146 
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N3393 
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16.8 
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1 
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0.784 








N2787 


17.2 


46 


0.39 


14.8 


6.7 
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0.30 


17.6 


21.3 


0.2 


1 
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N3384 
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59 


0.49 


14.5 


6.5 


1.9 


0.13 


/ 


/ 




1 


6.77 


0.479 


0.521 









Notes: Column (1), name of the galaxy. Column (2-4), the disc central surface brightness, scalelength, and ellipticity. Column (5-8), the bulge effective 
surface brightness, effective radius, Sersic index, and ellipticity. Column (9-11), the bar effective surface brightness, effective radius, and Sersic index. 
Column (12) and (13), apparent magnitude of the central source and the total galaxy. Column (14-17), luminosity fraction of the bulge, the disc, the bar 
and the central source. Luminosity parameters are in units of mag arcsec~ 2 , and scalelengths in arcsec. 
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Notes: Column (1), name of the galaxy. Column (2), Hubble type of the galaxy, "c" denotes the core elliptical galaxy, "p" denotes the pseudo-bulge, "p" 
denotes coexistence of the classical bulge and the pseudo-bulge. Column (3), Distance to the galaxy in units of Mpc. Column (4), mass and ltr error 
of the central black hole. Column (5), A"-band bulge luminosity in units of Lq k. Column (6) and (7), B-V and r-i color of the bulge. Column (8), 
bulge effective radius in units of kpc. Column (9), effective stellar velocity dispersion in units of km s , Column (10) and (11), the bulge stellar mass, 
calculated by K band M / L derived from B-V and r-i color. Column (12) and (13), the bulge dynamical mass calculated by the Sersic model and the 
isothermal model. All mass quantities are in units of M©. 
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Figure 2. 2MASS K" band image decomposition result. From left to right, we display the original image, model image, residual image, and surface brightness 
profiles of different components along the semimajor axis of the bulge, including the original image (black filled circles), the total model (grey thick line), the 
bulge (red dashed line), the disc (blue dot-dashed line), the bar (magenta dotted line) and the central compact source (green solid thin line). The thick scale bar 
represents 30 arcsec, except for NGC 224, which represents 300 arcsec. 
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Figure 2. 2MASS K band image decomposition result (continued). 
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Figure 2. 2MASS K band image decomposition result (continued). 



use the calibration by Bell et al. (2003): log (MJL K )=0A35(B-V)- 
0.356, and log (MJL K )=0349(r-i)-0336. We choose the SDSS r-i 
color because it is the most sensitive one to M/L (cf. Table 7 in 
in Bell et al. 2003). The value of extinction corrected B-V color is 
taken from the HYPERLEDA database. We use the total B-V color 
for elliptical galaxies and B- V color within the effective aperture for 
bulges as the approximation. In order to determine M/L of bulges 
more accurately, we also directly measure the r-i color within the 
bulge effective radius from the available SDSS image. For galax- 
ies with AGN component, the central 3 arcsec regions are removed 
in color measurement to avoid spectral contamination from AGNs. 
The formula of Bell et al. have statistical uncertainties of 0.1-0.2 
dex, thus the la error of logA/ s is adopted as 0.15 dex. 
The dynamical mass of bulges are estimated by 

M dyn = kR c a 2 JG, (4) 

where R a is the bulge effective radius, k is a model dependent di- 
mensionless constant, G is the gravitational constant. In the isother- 
mal model (<r» is a constant throughout the galaxy), we follow 
MH03 to use k = 3 instead of 8/3. In the more realistic Sersic 
model, k is a function of the Sersic index n, which is determined 
numerically. The details of the Sersic dynamical model is described 
in Appendix A. The la error of R c and a* are 10%-20% and 5%, 
the uncertainties of Sersic index An ~ 0.5 induce Ak/k ~15%. 
According to eq. (4), the la errors of log Ms in the isothermal and 
Sesic model are adopted as 0.1 dex (25%) and 0.15 dex (40%) re- 
spectively. 

We compare M s and Mi in Figure 3. M s calculated from the 
two colors are consistent very well, the systematic difference is 
much smaller than their errors. A/ d estimated by the Sersic model 
is systematic higher than Af d by the isothermal model, which is a 
natural result of the fact that k > 3 for all of our sample galaxies in 
the Sersic model (cf. Figure Al). MJM S is mass dependent, higher 
for masive bulges and lower for small bulges. The systematic dif- 
ference can be as high as 2.5 times. 

We also compare our results with the previous work in Figure 



4. Our L bl ,i.K and A/ d . iso are consistent with that given by MH03 for 
the elliptical galaxies, but are systematic (~0.3 dex) smaller than 
MH03 for the bulges. It may reflect the difference of decomposi- 
tion or data fitting methods used by BUDDA and GALFIT. However, 
the details and images of the two-dimensional decomposition in 
MH03 have never been published. As we have mentioned, the rib 
and R c are strongly coupled. The determination of M d , iso is sensitive 
to Tie, while L M-K is relatively robust for different decomposition 
programs. Our M 4Sh is similar with M d given by HR04, both are 
calculated under the assumption of Jeans equation, thus more reli- 
able than Af d , iso . The value of our M d-S er are consistent with HR04 
for the elliptical galaxies, proving the accuracy of our method, but 
are systematic (~0.3 dex) smaller than HR04 for the bulges. The 
difference are due to our decomposition for the SO and some ellip- 
tical galaxies, while HR04 treat them as a whole, thus overestimate 
the bulge mass. The 0.3 dex deviation indicates about half of the 
luminosities of SO galaxies come from disk component, consistent 
with our decomposition results (cf. column 14 and 15 in Table 1). 



3 RESULTS 

As in H08, we use two kinds of bisector linear fitting methods, 
"X 2 " (Press et al. 1992) and "AB" (Akritas & Bershady 1996), to 
determine the black hole-bulge relations in the form of log M bb = 
a + P(logx — xq), where a and @ are the intercept and slope of 
the relation, x is the value of the bulge property concerned, xo is a 
chosen constant. 

The fitting results of Mim-Mm correlations are presented in 
Table 3. eo is the intrinsic scatter of the correlation given by the \ 2 
method. The parameters of relations derived by two fitting methods 
are almost the same. 

The M bh -L bulK relations are shown in Figure 5. The classical 
bulges and pseudo-bulges have distinct correlations. In the residual 
diagram of the Af bh -L bll i.K relation, we discriminate the core ellipti- 
cal, normal elliptical, classical bulges, and pseudo-bulges. There is 
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Figure 3. Comparison of A/ S and M&, The solid lines denote the linear fit. The errorbars are not shown. 




no obvious type dependence in the residuals. We also try to fit the 
M bb -L M _ K (the total _ftT-band luminosity of the galaxy) relation for 
the classical bulges, find a looser relation (Figure 6), which confirm 
the previous results that bulges are better tracer of SMBHs than the 
whole galaxies. 

The Mbh-Af, and M bh -A/ d relations are shown in Figure 7 and 
8. In all the cases, pseudo-bulges follow independent relations with 
similar slope and over 1 dex smaller intercept comparing with the 
relations for classical bulges. The number of pseudo-bulges in Af bb - 
M s _„ diagram is only three, we do not fit them but note they also 
locate below the classical bulges. The two galaxies harbor both the 
classical and pseudo-bulges just locate between the relations of two 
type of bulges, reflect the mixed nature of their bulge properties. 
The different black hole-bulge relations obeyed by the two type of 
bulges are emphasized in Figure 9. The statistical significance of 
the difference are over 3<r limit. 

We also fit these black hole-bulge relations for core elliptical 
galaxies and compare the results with that of the other classical 
bulges (Figure 10). Their M bb -L bu]x and Af bh -Af s relations seem to 
be slightly shallower than the others, while their M bh -M d relations 
are consistent within la limit. 

Our results for classical bulges are compared with various past 
work (Table 4). The slope j3 ~ 1 are consistent with the previous 
value, but the intercept a are 0.1-0.3 dex larger. These systematic 
differences are partly due to the more accurate decomposition for 
bulges, and partly due to the removal of pseudo-bulges in fitting. 



Although our sample is the biggest one, the intrinsic scatter of our 
Afbh-Lbui.K relation is still larger than the M bh -a* relation, while our 
Af Ml -M d relations are as tight as the Af bh -cr, relation. M d is a better 
tracer of A/ bh than L w > it is as good as a* but more robust for core 
elliptical galaxies. 



4 CONCLUSIONS AND DISCUSSION 

We have examined the black hole-bulge correlations in a sample of 
58 galaxies, especially the bulge type dependence of the correla- 
tions. The main results are as follows: 

(1) The improved versions of M bh -M bul relations for elliptical 
galaxies and classical bulges are: 

log M bb = (8.38 ± 0.05) + (0.97 ± 0.08) log(L bul . K /lO ia9 L ,K) 
= (8.24 ± 0.05) + (1.07 ± 0.09) log(M slc /10 10 ' 5 M Q ) 
= (8.46 ± 0.05) + (0.90 ± 0.06) log(M dyn /lO 11 M ), 

with intrinsic scatter eo =0.36, 0.32, 0.27 respectively. 

(2) The pseudo-bulges follow different M bh -M blll relations: 

log Af bb = (6.48 ± 0.14) + (1.15 ± 0.37) log(JW/lQ lo - 3 L 0|K ) 
= (6.77 ± 0.13) + (1.36 ± 0.40) log(A/ slc /10 10 ' 3 M Q ) 
= (6.67 ± 0.10) + (0.85 ± 0.19) log(Af dy „/lO ia3 M ), 

with intrinsic scatter eo =0.36, 0.21, 0.19 respectively. 
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Figure 5. Top: The Af bb -L bll i, K relations. Bottom: The Residuals between 
the observed M b i, and the predicted value from the best-fit Mbh-LwjE rela- 
tion of elliptical galaxies and classical bulges. Thick solid and dashed lines 
denote best-fit and la uncertainties of the relation for classical bulges, the 
thin lines denote the corresponding relation for pseudo-bulges. The dash- 
dotted line denote the best-fit result given by MH03. The filled circles de- 
note the core elliptical galaxies, the stars denote the other elliptical galax- 
ies, the open squares denote the classical bulges, the open circles denote the 
pseudo-bulges, triangles denote the galaxies harbor both the classical and 
pseudo-bulges, the points with dashed errorbars denote the edge-on galax- 
ies. 
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Figure 6. The M bb -L, olK correlation. The symbols denote as the same as in 
Figure 5. 



Table 4. Comparison with previous results [log M b h= a + /3(log x — xq)] 



X 


Xo 


Q 


P 




Nt 


Ref. 


-^ b ul.K 


10.9 


8.21±0.07 


1.13±0.12 


0.31 


27 


MH03 


-k b ul,K 


10.91 


8.29±0.08 


0.93±0.10 


0.33 


22 


G07 


-k b ul,K 


10.9 


8.38±0.05 


0.97±0.08 


0.36 


47 


this work 


-^ b ul,V 


10.33 


8.41±0.11 


1.40±0.17 




41 


L07a 


-^ b ul,V 


11.0 


8.95±0.11 


1.11±0.18 


0.38 


38 


G09b 


M, 


10.5 


8.24±0.05 


1.07±0.09 


0.32 


47 


this work 


M d 


11.0 


8.20±0.10 


1.12±0.06 


0.3 


30 


HR04 


M d , s „ 


11.0 


8.46±0.05 


0.90±0.06 


0.27 


47 


this work 


M iim 


10.9 


8.28±0.06 


0.96±0.07 


0.25 


27 


MH03 


M d , iso 


10.9 


8.61±0.05 


0.88±0.06 


0.27 


47 


this work 




2.30 


8.28±0.05 


4.06±0.28 


0.27 


39 


H08 




2.30 


8.23±0.08 


3.96±0.42 


0.31 


25 


G09b 



< the number of the sample galaxies for fitting. 



At a fixed bulge mass, Mm, in pseudo-bugles are on average 
over one magnitude smaller than that in classical bulges. 

(3) Bulge dynamical mass M d is a better tracer of A-/ bb than 
L bu i and M s . The Af bh -A/ d relation is as tight as the Af bb -<T* relation, 
both have intrinsic scatter eo < 0.3 dex. 

(4) The core elliptical galaxies obey the same M^-Ms relation 
with the classical bulges. 

We have invalidated the assumption used by Gadotti & Kauff- 
mann (2009), i.e. the pseudo-bulge do not follow the same Af bb - 
M s relation with the classical bulges. For a pseudo-bulge with a 
give M s , the corresponding Af bh is probably much smaller than they 
have assumed, making their deviation in the A/ bb -a» diagram even 
larger (cf. Figure 1 in their paper). Therefore, their results will be 
strengthened rather than weakened, and confirm the conclusion of 
H08 again. 

Our work have confirmed the results of Greene et al. (2008), 
that the psudo-bulges have much smaller A/ bh than in the classical 
bulges. Their sample pseudo-bulges are distant and growing. The 
growth of SMBHs in pseudo-bugles seems not as efficient as in 
the classical bulges. The fueling gas for SMBHs accretion may be 



short in supply, or be consumed by competitive mechanism, e.g., 
formation of nuclear star clusters (Nayakshin et al. 2009), or their 
growth timescale is longer than the Hubble timescale. 

The mass of giant elliptical galaxies can be increased by 25% 
since z ~ 1 by dry stellar mergers (e.g., Naab et al. 2007). Sup- 
pose the SMBHs coalescent finally, and no significant star forma- 
tion or accretion take place in dry mergers, the M bb /M bll i should 
preserve since high redshift. If we believe the core elliptical galax- 
ies are the product of dry mergers, their slightly shallower slope in 
the M bh -M h relation indicate that cosmic evolution of Af bb /M s may 
exist. However the average Af bb /Af d ratio is the same, this differ- 
ence may related to the mass contribution from the dark matter, or 
to the selection bias in the high mass end of correlations (Lauer et 
al. 2007b). We need larger sample and more reliable Af bb measure- 
ments to check this suggestion. 

Our new black hole-bulge relations have important implica- 
tions for the studies of SMBH demography and coevolution mod- 
els of black holes and galaxies. We will explore these problems in 
forthcoming papers. 
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Table 3. Black hole - bulge correlations [log Mt,h= ct + /3(log x — xg)] 



X 


xo 


a x 2 


/3 X 2 


eo 




/3ab 








classical bu 


ges 








10.9 


8.38±0.05 


0.97±0.08 


0.36 


8.38±0.05 


0.97±0.09 


L to t, K 


10.9 


8.16±0.07 


1.06±0.13 


0.45 


8.16±0.07 


1.06±0.12 




10.5 


8.24±0.05 


1.07±0.09 


0.32 


8.24±0.05 


1.08±0.11 




10.5 


8.25±0.06 


1.17±0.14 


0.29 


8.25±0.06 


1.18±0.16 


Md.Se, 


11.0 


8.46±0.05 


0.90±0.06 


0.27 


8.46±0.05 


0.90±0.07 


Md. is „ 


10.9 


8.61 ±0.05 


0.88±0.06 


0.27 


8.61±0.05 


0.88±0.06 








pseudo-bul 


ges 








10.3 


6.48±0.14 


1.15±0.37 


0.30 


f. AS4-() no 


1 .27±0.34 




10.3 


6.77±0.13 


1.36±0.40 


0.21 


6.79±0.11 


1.49±0.46 


M d , s „ 


10.3 


6.67±0.10 


0.85±0.19 


0.19 


6.68±0.08 


0.89±0.19 


M d . iso 


10.3 


7.01±0.13 


0.81±0.19 


0.22 


7.01±0.11 


0.82±0.16 






core elliptical g 


alaxies 






-^bul.K 


10.9 


8.59±0.12 


0.73±0.21 


0.38 


8.59±0.10 


().73±0.19 




10.5 


8.44±0.15 


0.81±0.22 


0.36 


8.43±0.14 


0.82±0.25 


M d , s „ 


11.0 


8.48±0.10 


0.86±0.17 


0.29 


8.47±0.08 


0.87±0.16 


M d . iso 


10.9 


8.60±0.09 


0.88±0.16 


0.29 


8.60±0.07 


0.88±0.16 
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'"SC^dyn sun' [ isotherma| ] 

Figure 8. The M\trMt relations. The dash-dotted line denote the best-fit 
results given by HR04 (top) and MH03 (bottom), other symbols denote as 
the same as in Figure 5. 



APPENDIX A: DYNAMICAL MASS ESTIMATION IN THE 
SERSIC MODEL 

The Sersic dynamical model is described as follows. 

(1) The galaxy surface brightness distribution 1(f) is a Sersic 
profile, 



I(r) = I eM-C n (r/r c ) 1/n ], 



(Al) 



where r is the Sersic index, r c is the effective radius, C n — 2n ~ 
0.324 is a constant. 

(2) The mass-to-light ratio (M/L) is a constant throughout the 
galaxy, 



M 

T 



J(r) = / P (f) 



2r'dr' 



(A2) 



whrer p(r) is the stellar density. 

(3) The gravitational potential <&(r) is dominated by stars, the 
other mass components (e.g., gas and dark matter) are negligible. 



d$(r) _ GM(r) 
dr r 2 



7.5 

7- 
6.5 
* 6 
5.5- 

5- 

4.5 % 

4 1 1 1 1 1 

1 2 3 4 5 

n 

Figure Al. The calculation result of k(n) in the Sersic model. The solid 
line is for <r e , the dashed line is for a c . 

where M(r) = f Q r 4irr' 2 p(r')dr' is the enclosed mass within ra- 
dius r. M(po) is the total mass. 

(4) The stellar velocity dispersion o(r) is isotropic, follow the 
Jeans equation, 

d[p(r)q 2 (r)] 



-p(r[ 



d$(r) 



(A4) 



dr dr 

Given M, r e , and n, we can derive <r(r) by solving eq. (Al- 
A4). The projected flux-weighted stellar velocity dispersion a p at 
radius r is 

a f x o 2 (r')p(r')(r' 2 - r^-^^r'dr' 

(M/L)Hr) ■ (A5) 

The observed flux-weighted central stellar velocity dispersion 
within a slit aperture of length 2ro is 

(A6) 



al{ro) _I-I(r)4(rW 



Jo" /(r)dr' ' 

The observed flux-weighted central stellar velocity dispersion 
within a circular aperture of radius rn is 

/;° r'I(r)a 2 (r')dr' 



o c (r ) 



Io° r'I(r)dr' 



(A7) 



a, can be defined as a* = <T e (r e ) or cr* = o" c (r e /8), i.e. the 
so-called "effective velocity dispersion" or "central velocity dis- 
persion". They are roughly equivalent with minor difference (cf. 
discussion in H08). We calculate k(n) in eq. (4) for both cr e and 
<t c , the results are shown in Figure Al. In this paper, we use the 
definition of cr e . 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 
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